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Abstract. We study inhomogeneous semilinear parabolic equations with 
source term / independent of time u, - Au + + f(x) on a metric mea- 
sure space, subject to the conditions that f{x) > and m(0, x) = (f{x) > 0. 
By estabHshing Harnack-type inequaHties in time f and some powerful 
estimates, we give sufficient conditions for non-existence, local exis- 
tence, and global existence of weak solutions. This paper generalizes 
previous results on Euclidean spaces to general metric measure spaces. 
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1. Introduction 

In recent years, the study of PDEs on self-similar fractals has attracted 
increasing interest, see for example [|7] |8] |9l [131 [II]- In this paper we in- 
vestigate a class of nonlinear diffusions with source terms on general metric 
measure spaces. Diffusion is of fundamental importance in many areas of 
physics, chemistry, and biology. Applications of diffusion include: (1) Sin- 
tering, i.e. making solid materials from powder (powder metallurgy, pro- 
duction of ceramics); (2) Catalyst design in the chemical industry; (3) Steel 
can be diffused (e.g. with carbon or nitrogen) to modify its properties; (4) 
Doping during production of semiconductors; (5) The well-known Black- 
Scholes Model in Financial Mathematics that is closely related to option 
pricing can be transformed to a parabolic equation. 
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Let (M, d,iu) be a metric measure space, that is, (M, d) is a locally com- 
pact separable metric space and yti is a Radon measure on M with full sup- 
port. We consider the following nonlinear diffusion equation with a source 
term / on {M,d,jj): 

Ut = Au + + fix), t > and x e M, (1.1) 

with initial value 

u{0,x) = (p{x), (1.2) 
where p > I and /,(/?: M —> R are non-negative measurable functions. 
With an appropriate interpretation of weak solutions of (11.11) on {M,d,jj), 
we shall investigate the non-existence (or blow-up) of solutions, the local 
and global existence of weak solutions to (ll.ll) - (ll.2l) . as well as the regu- 
larity of these solutions. Although we were partially motivated by a series 
of earlier papers [H] |9] HO] [H] [151 [HI [HJ, there are new ideas in this paper. 
In particular, we have used the theory of heat kernels on metric measure 
spaces. 

Recall the definition of the heat kernel which will be central to our ap- 
proach. A function k{-, ■ , •) : M x M ^ called a heat kernel if the 
following conditions (fcl)-(M) are fulfilled: for //-almost all (jc, j) e MxM 
and for all t,s > 0, 

(kl) Markov property: k(t,x,y) > 0, and j^k(t,x,y)dfx(y) < 1; 
(kl) symmetry: k(t,x,y) = k(t,y,xy, 

(k3) semigroup property: k(s + t,x,z) = k(s, x, y)k(t, y, z)diJL{y); 
(M) normalization: for all / e L^{M,ju) 



lim I k(t, X, y)f{y)dix{y) = f(x) in the L^(M, //)-norm. 

We assume that the heat kernel k(t, x, y) considered in this paper is jointly 
continuous in x, y, and hence the above formulae in (fcl)-(M) hold for every 
(x,y) e MxM. 

Two typical examples of heat kernels in R" are the Gauss-Weierstrass and 
the Cauchy-Poisson kernels: 

^(^'-^'3^) = 7^-^777 exp 



{Antfl^ \ At 



k{t,x,y) = — |1 + — 



t" \ t^ 



-^"^'^'^ < r(kn + i)) 



c = 



^(n+l)/2 ^ 



Jointly continuous sub-Gaussian heat kernels exist on many basic fractals, 
for example, on the Sierpmski gasket, see Barlow and Perkins [5], and on 
Sierpmski carpets, see Barlow and Bass [|3l |2l. For other fractals see [[T2l 
[T3l . For non-sub-Gaussian heat kernels, see [141161. 
A heat kemel k is called conservative if it satisfies 

(^5) conservative property: J^k(t,x,y)diJ.(y) = 1, for all f > and all 
X e M. 

We will also assume that the heat kemel satisfies the following estimates 
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(^6) two-sided bounds: there exist constants a,/3 > such that for all 
t > and all x,y e M, 

1 ^ ( dix,y) \ ^ 1 ^ ld(x,y)\ 

where Oi and O2 are strictly positive and non-increasing functions 
on [0, 00). 

It turns out that the parameter a in (11.31) is the fractal dimension, and /3 is 
the walk dimenison of M, see [11]. 

Two-sided estimates (11.31) hold on various fractals where 

0,(^) = C- exp(-c;//^-^^) (for all s > 0) 

for constants C,, c, > (z = 1,2) and /3 > 2 is the walk dimension. 

To prove the regularity of solutions, we need to assume that the heat 
kernel k is Holder continuous in the space variables: 

(kl) Holder continuity: there exist constants L > 0, v > 1 and < cr < 1 
such that 

\k(t,xi,y) - k(t,X2,y)\< Lt~''d(xuX2T 
for all? > and all xi,X2,y e M. 
Given a heat kernel k, the operator A in (11.11) is interpreted as the infini- 
tesimal generator of the heat semigroup {^r},>o in := L}{M,ij). Thus we 
let 

K,g{x)= r k{t,x,y)g{y)dix{y){t>Q,geL\ (1.4) 
Jm 

and define A by 

^g = lim (in L^-norm). (1.5) 

ao t 

Observe that {^f}f>o is a strongly continuous and contractive semigroup in 
L^, that is, for all s,t >0 and all g e L^, 

Ks^t = K.Kr, (1.6) 
lim||i^,^-^||2 = 0, 

WKM, < \m, (foralll <^<oo). 

The domain of the operator A is dense in L^. 

A function u(t, x) is termed a weak solution to (ll.ll) - (ll.2l) if it satisfies the 
following integral equation 

u(t, x) = Kt(p(x) -\- I Krf{x)dT+ I K,.rU''{T,x)dT, (1.7) 
Jo Jo 

where K, is the heat semigroup defined in (11.41 ). 

The structure of this paper is as follows. In Section |2l we show the non- 
existence of weak solution to (ll.ll) - (ll.2l) . In Section [31 we obtain sufficient 
conditions for the local and global existence of solutions for a range of 
parameters p, source terms / and intial values (p. The critical exponents p 
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depend only on the fractal dimension a and the walk dimension /3. Finally, 
in Section m we investigate the Holder continuity of weak solutions. 

Notation. The letters C, C,(/ = 1,2,...) denote positive constants whose 
values are unimportant and may differ at different occurrences. 



In this section we give sufficient conditions for the non-existence of es- 
sentially bounded solutions. The exponents p = I + /3/a (where a,/3 > 0), 
and p = a/{a - /3), where a > /3 > occur in the heat kernel bounds (11.31 ). 
play a crucial role in our analysis, see Theorem 12.21 First, we establish 
Lemma [Zn where condition (k6) is our only assumption on the heat kernel 
k (we do not need the conservative property of k at this stage.) 

The following properties the functions Oi and O2 in condition (^6) may 
or may not hold: there exist positive constants a,, bi and such that, for all 
s,t>0. 



Note that if (|2.1I) holds, then < ai < 1 by letting s = and using the fact 
that 02(0) > Oi(0). Without loss of generality, we may assume that ^2 > 1 
in (12.11) . since if (12.11) holds for some 02 < 1, it also holds for any constant 
^2 > 1 by the monotonicity of O2. 
The Gauss-type functions Oi and O2 



for constants y > and C, > 0(1 < / < 4) satisfy properties (|2.11 )- (|2.31 ). 
The Cauchy-type functions 



for constants 7 > and C, > 0(z = 1,2), satisfy properties (12.11 ) and (12.21 ). 
but not (1231) if p > I. 

Condition {k6) and inequality (12.11) lead to the following key lemma. 

Lemma 2.1. Assume that the heat kernel k satisfies condition (k6) and ( I2.il) . 
Then, for all non-negative measurable functions g on M and for all t > 
0,x e M, 



2. NON- EXISTENCE OF SOLUTIONS 



Oi(5) > ai02(a25), 
^lis + t) > bi^2(b2S)^2(b3t), 
0{'(5) > Ci02(C25). 



(2.1) 
(2.2) 
(2.3) 



Oi(5) = Ciexp(-C25^), 
02(5) = C3 exp(-C45^), s > 0, 



(2.4) 



Oi(^) = Cl(l+5)-^ 
02(5) = C2 (1 + s)-^ , s>0 



(2.5) 



Krg(x) > AiKBrg(x), 



(2.6) 




(2.7) 
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where Ai = aia^" < 1,^2 = aia^"{\ — a^~) < 1 and B = < 1. 
Consequently, for all non-negative measurable functions (p, 

KMX) + Krg(x)dT > A [KsM^) + tKj,^,g(x)] , (2.8) 

where A = min jA^ A2} < 1 and Bi = = a'^^ . 
Proof. It follows from condition (^6) and (12.11) that 

Ktg{x) = I k{t,x,y)g(y)dfi(y) 
Jm 

which gives that, using (^6) again, 

Ktgix) > aia2" \ k{a~^t,x,y)g{y)dn{y) 
= aia2°'K^H^fg(x) = AiKBtg(x), 

proving (|2.6I ). 



To show (|2.7I) . we see from (12.91 ) that for all t e t], using the mono- 
tonicity of O2 and condition (^6), 

r 1 / d(x,y)\ 
Krg(x) > a I J^^(t>2\a2 ^^ip \g{y)dix{y) 



> 



r 1 ^ f d{x,y) \ 



g{y)dix{y) 



a2^" I k{a2'^t,x,y)g{y)dix{y) 
aia2^"KB2tg(x). 



> ai 

-2a 



Therefore, 



Jo 







> 


f 




Ja-/t 






> 


r 




Ja-/t 







aia2^"KB2fg{x)dT 



proving (12.71 ). 

Finally, replacing t by Bt, we see from (12.61) that Kst^pix) > AiK^itcpix) = 
AiKb^,(P(x), and thus 

Kt<p(x) > A.KbMx) > AKsMx)- (2.10) 
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Adding (O and (llJOl) . we obtain ■ 

Lemma 12.11 gives the following estimate (12.111) that plays an important 
role in proving the non-existence of global bounded solutions. 

Theorem 2.2. Assume that the heat kernel k satisfies conditions {k6) and 
l\2.1\) . Let u(t, x) be a non-negative essentially bounded solution of U.7^ in 
(0, T) X M. Then, for all (t, x) 6 (0, T) x M, 



where Bi = a^ as before, and Ci depends only on p (and in particular is 
independent ofT, <p and f). 

Proof. Observe that by condition (k\) and using a weighted Holder 
inequality, for all t > 0,x e M and for all non-negative functions g. 



t'/^P-'^KsMx) + tP'^P-'^KB.Jix) < C 



(2.11) 




It follows from (O) and (l2J0l) that 




(2.12) 



From (|1.7I) and (|2.81 ). we see that 




> A{Kj,Mx) + tKB,J{x)). 



(2.13) 



Starting from (12.131 ). we shall apply (12.121 ) repeatedly to deduce the desired 
inequality (|2.1 11) . Indeed, we obtain from (12.121) and (12.131) that, using the 
semigroup property (II.6I 1 of {K,}f^Q and the elementary inequality (a + by > 
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aP + bP for all j9 > 1 and a,b>0, 

u(t,x) > A I [KBi(,-r)UiT, x)Y dr 
Jo 

> A f [Ks,^,.r){A(Ks,r<P + TKs,rf)}{x)YdT 

Jo 

= AP^' [ [KBM^) + rKB,tf{x)rdT 
Jo 

= AP^' UkbM^))' + -T^t'^^KBjm 

I 1 +p 

Repeating the above procedure, we obtain that for all n > 1 , 



+ 



(1 + p)P"'\\ +P + p^)P"'' •■• (1 +;? + •■• + 

t'^p^-^p^KBtmy" 



(1 + p)P {I +P + p^)P • • • (1 + p + • • • + p")^ 
It follows that 

n 

<u(t,x)P'" Y]il+p + --- + p'-Y\ (2.14) 

i=2 

' - l)/(p- l)p" - l)/(p- 1 )p" KBtf(x) 

n 

<u(t,x)P " ]~[(1 +;? + ■•■+ '. (2.15) 

!=1 

Since 



iog]~[(i + /? + ••• + /-!)"' < 2 - logoy') < +00, 

i=2 (=2 

logf](l +;? + •■■+/)"' < _2-log((/+l)y)<+c«, 

(=1 !=1 ^ 

and that u{t, x) is essentially bounded on (0, T) x M, we pass to the limit as 
n ^ oo in (12.141) and (|2.15l) . and conclude that 

t"^P-'^KBtip(x) < Ci/2, (2.16) 
tP/(p-^)KBtf(x) < Ci/2, (2.17) 

for some Ci > 0. Adding (12.161) and (12.171) . we obtain (12.111) . ■ 
We are now in a position to obtain the main results of this section. 
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Theorem 2.3. Assume that the heat kernel k satisfies conditions (k6) and 
( \2.1\l . Then the problem di.iD -f [772]) does not have any essentially bounded 
global solution in each of the following cases: 

(i) z/p < 1 + ^ and if either ip{x) ^ or f{x) ^ 0. 

(ii) ifa</3 and if f{x) ^ 0; 

(iii) if a > /3 and p < ^(> 1 + f ) and if f{x) ^ 0. 

Proof. We prove the results by contradiction. Assume that u{t, x) is a 
non-negative essentially bounded global solution. Replacing Bit by t, we 
see from (12.111) that for all x e M and t > 0, 

t"^P-'^Ktip{x) + t"'^P-'^K,fix) < Ci, (2.18) 

where < Ci < oo is independent of (p and /. 

Proof of Case (i): If (p{x) ^ 0, we see from (k6), using Fatou's lemma, 
that 

liminf t"'^KM^) > liminf r 01(^^)^(3;)^^^) 

> C2. 

where C2 = 1 if ||^||i = 00, and C2 = Oi(0)||^||i if \\(p\\i < 00. However, as 
1 /(p - 1) > a//3, this is impossible by using (12.181) . Hence, (fLT])-(fr2l) does 
not have any global essentially bounded solution. 



If f{x) ^ 0, observe that u{t + to, x) is a weak solution of (11.71) with initial 
data (fix) = u(t[),x). We may find to > such that u(to,x) ^ 0. Repeating 
the above argument, we again see that (|1.11 )- (|1.21 ) does not have any global 
essentially bounded solution. 

Proof of Case (ii): Observe that by (11.71 ) and (12.71) . 



uit,x)> [ KJ{x)dT > A^tK^Jix). (2.19) 
Jo 

We distinguish two cases: or < /3 and a = /3. 

• The case a < /3 . It follows from (12.191) and (^6) that 

liminf t^""^^-^u{t,x) > A2liminf t""^KB„f{x) 

> C3, (2.20) 

where C3 = 1 if ||/||i = 00 and C3 = A25^''^^Oi(0) > if ||/||i < 00. 
However, since u is globally essentially bounded and a//3 < 1, we see 

liminf /"'^^-^uit, x) = 0, 

a contradiction. 
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• The case a = fi. For t > 1, it follows from (^6), ( 12.11) and the mono- 
tonicity of O2 that 

u(t,x) > J Krf{x)dr 



> a, I dr I T-'(^2{a2^^^\f(y)diu(y) 



d(x,y) 



fdr f 

a, [ T~'dT r (:>2{a2d{x,y))f{y)dii{y). (2.21) 
Ji Jm 



Since f{x) ^ 0, we can find a point x e M such that 

^2{a2di,x,y))f{y)dn{y) > 0. 



M 



Passing to the limit as ? ^ 00 in (12.211) this contradicts that u is globally 
essentially bounded. 

Proof of Case (Hi): It follows from (12.181) and (^6) that 

liminfCi/'''/^'P/^''-i^ > \immi f'l^K J {x) 



> C4, 

where C4 = 1 if ||/||i = 00 and C4 = <I>i(0)||/||i if ||/||i < 00. However, this 
is impossible since ^ - < 0. The proof is complete. ■ 

In Theorem [231 we do not know in general if there exists any essentially 
bounded global solution for two the critical cases p = I + - {a,P > 0) and 
p = ^^{a>l5>0). 

However, Theorem 12.3 1 (/) may be improved to include the critical expo- 
nent p = 1 + ^ under further assumptions (12.21 ) and (12. 3t on the heat kernel 
k. We first need the following property. 

Proposition 2.4. If<^2 satisfies A2.2\) . then for all t > and all x,y e M, 

(!)2(dix,y)t-''^) 



> 



b,(^2(b3d(y,0)t-"f'), (2.22) 



O2 (M(^,0)ri//^) 
where the constants bi (i = 1, 2, 3) are as in l\2.2\l . 

Proof. Since O2 is strictly positive and decreasing on [0, 00) and d(x, y) < 
d{x, 0) + d(y, 0), we have 

(^2(d(x,y)t-^'^) > (^2(d(x,0)t-^'^ + d(y,0)t-^''^) . (2.23) 



It follows from (|Z21) that 
<D2 (dix, 0)t-^'^ + d(y, Or^^^) > bi(^2 (b2d{x, QX^"^) O2 (b^d{y, 0)r'/^) , 
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which combines with (12.231) to give (12.221) . ■ 

Theorem 2.5. Assume that the heat kernel k satisfies conditions (k5), (k6) 
and d2.iD . ( 12.21) and ( 12.31) . Then ( Ii.il) -( I772]) does not have any essentially 
bounded global solutions if p < I + ^ and if either <p{x) ^0 or f{x) ^ 0. 

Proof. In view of Theorem 12.31 (z) it is enough to consider the critical 
exponent p = I +13 1 a. We only consider the case ip{x) ^ (the case f{x) ^ 
may be treated in a similar way). Then (12.181 ) becomes 

t"'^K,(p(x) + t^^""^KJ(x) < Ci. 

From condition {k6) 

f ip{y)dix{y)<C2, (2.24) 
Jm 

where C2 = Ci/Oi(0). For any to > 0, the function v{t,x) = u{t + to,x) 
is a weak solution to (11.71) with initial data ip(x) = u(to, x). Repeating the 
procedure of (12.241 ). we have that for all to > 0, 

f u(to,y)dij{y) < C2. (2.25) 
Jm 

We claim that there exist positive constants 7, p possibly depending on to 
and (f such that, for all x e M, 

u{tQ, x) > pk{y, X, 0). (2.26) 

To see this, observe that 

^2{d{x,0)y-'l^)>k{y,x,0)y"ll', 
and thus, using (12.221) and setting 7 = (aiZ?2)~^?(), 

O2 [a2d{x, y)t-^"^) > ^i02 (a.b^diy, 0)?"'^^) O2 (aiM(x, 0)?"^^^) 
> br^2 (aMy, 0)t-"^) k{y, x, 0)y"l''. 

Using (O) and (IXTT) . 

u{tQ,x) > I k(to,x,y)(p(y)diu(y) 
Jm 

> t-"'" [ ^,(d(x,y)t-"'')<f(y)diu(y) 

Jm 

> a.t-"'" f (!>2(a2dix,y)t~'"')<p(y)diu(y) 

Jm 

^1 k{y,x,0) J (^2[alb^d{y,0)to"')viy)d^l(y), 
and hence, inequality (12.261) holds by setting 

^1 O2 (a.b^diy, 0)t-"'') ^(y)dM(y), 

proving our claim. 
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Consider v{t, x) = u(t + to, x) such that u(to, x) ^ 0. Applying (12.261) . we 
obtain 

v(t,x) > I k(t,x,y)u(t(i,y)diu(y) 
Jm 

> p I k{t,x,y)k{y,y,0)dn{y) 

JM 

= pk(t + y,x,0). 



which yields that, using (11.71) . (kS) and Fubini's theorem, 

Jv(t, x)djjL{x) > I dix{x) I dT I k{t - r, x, y)v''{T, y)diJ.(y) 
M Jm Jo Jm 

v'\T,y)diu(y) 



fdr f 

Jo Jm 



IM 

> pP [ dr [ kPiT + y,y,0)dfi(y). (2.27) 
Jo Jm 



As ;? = 1 + /3/a, we see from (12.31) and (^6) that 

F(T + r,y,0) > (T + r)-^'+"'^^^o^(j();,0)(T + 7)-'/^) 

> Ci(T + y)-(i+''/^^02 (c2d(y, 0)(T + y)-'''') 

= c,c-^"{T + yr'[c-/iT + y)r'^^2 (cidiy, 0)(t + y)-'"") 



> c,c-''{T + yr'k(c-/ir + y),y,0), 



'2 V' ' "■V^2 

which combines with (|2.27l) to give that 



r v{t,x)dn{x)>CiC-^"pP r (t + 7)-'Jt. (2.28) 
Jm Jo 



Passing to the limit as f ^ oo, we conclude that 

v{t,x)diJL{x) — > oo, 



/ 

Jm 



IM 

which contradicts (12.251) . ■ 

We note that our results agree with the earlier ones where M = W and 
is Lebesgue measure, and where the heat kernel k is the Gauss-Weierstrass 
function (so that A is the usual Laplacian), see IfTUlfTSlfTO fTTl. See also Q 
where M is a fractal and jj. is or-dimenisonal Hausdorff measure, and where 
k is the Gauss -type heat kemel on M. 



3. Existence OF SOLUTIONS 



In this section we give sufficient conditions for local existence and global 
existence of weak solutions. 
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Theorem 3.1 (Local-existence). Suppose that the heat kernel k satisfies 
(k6). Let b{t) be a continuously differentiable function on [0, Tq) satisfying 



b'{t) = bP{t) 
with initial value b{0) = I- If 



■ r WKrf 

Jo bir 



f\\o 
) 



■dT + \\KMl 



Jo Jo bir) 



■dT + \\K,. 



ds < 



(3.1) 



(3.2) 



then di.il) - 4772]) has a non-negative local solution u e L°°{(0, T), M)for all 
< T < To, provided that \\(p\\oo < oo. 



Remark 3.2. By Peano 's theorem, there exists some Tq > and some con- 
tinuous differentiable function bit) such that (123) holds in [0, Tq). Clearly, 
such a bit) is non-decreasing in [0, Tq). On the other hand, condition d3.2|] 
may be verified for some specific cases. For example, if ikS) holds and if 
f = 0,(p = C > 0, then 

bit) = [l-ip-l)CP-'t]-"^''-'^ 

satisfies ( \3.1\i in [0, Tq) where Tq = ip - 1)~^C"^^"'\ and di.2D also holds. 
As an another example, let f = 1,^ = and assume ik5) holds. Then, for 
p = 2, we see that bit) = 1 / cos t satisfies di.iD for t e [0, n/2), and that 
d3.2D holds. 



Proof. Define 



ait) = bit) 



' \\Krf\\c 

bir) 



■dT. 



Note that a(0) = and ait) > for ? e [0, Tq). Incorporating this into (13.11 ). 
we get 

b'it) = bit)[ait) + bit)\\KMLY-\ 

Moreover, 

b'it)ait) 



a'it) = \\K,f\U + 



bit) 



= \\KJ\\^ + ait)[ait) + bit)\\KMU 



Together with the initial conditions, these differential equations are equiva- 
lent to 

ait) = f WKrfWoodT-^ f aiT)iaiT)-^biT)\\Kr(p\W'dT, (3.3) 
Jo Jo 

bit) = 1+ r biT)iaiT)-^biT)\\Krifi\U"-'dT. (3.4) 
Jo 

Let 'K be the family of continuous functions u satisfying 

K,(pix) < uit, x) < ait) + bit)K,(pix) for all it, x) e [0, Tq) x M. (3.5) 



Jo Jo 



x)dT. 



INHOMOGENEOUS PARABOLIC EQUATIONS 

Define 

Tu{t, x) = Ktip(x) + 

We claim that ifue'H, then Tu e "7/, that is, 

Kt(p(x) < Tu{t, x) < a(t) + b{t)Ktip{x) iO<t<To,xeM) 
Observe that, using (^1), 

r K,.r[a(T) + b{T)Kripnx)dT 

Jo 

= [ dT [ kit-r,x,y)[air) + b(T)Kr<p(y)Ydniy) 
Jo Jm 

< [ [aiT) + bir)\\Kr^\\^r-'[a(T) + biT)KMx)]dT. 
Jo 

It follows from and (1331) that 

ru(t,x) < Ktip{x)+ r Krfix)dT+ f Kt-r[a{T) + biT)Kr<pf(x)dT 

Jo Jo 

< \f \\Krf\LdT+ [ a{T)[a(T) + b{T)\\KMUP-'dT 

Jo Jo 
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(3.6) 
(3.7) 



+ 



1+ r b(TmT)+b{T)\\Kr<f\UP-'dT 

Jo 



KMX) 



= a{t) + b{t)Ktip{x) 

using (13.31) and (13.41) . so (13.71) holds, proving our claim. 
For n = 0, 1, 2, • ■ • , define 

uoit, x) = Kt(p(x), 

Un+l(.t,x) = Tu„(t,x). 



Using (13.61) inductively, it follows that the sequence {«„(?, x)} is non-decreasing 
in n, and, for all n > and all x e M,t e [0, Tq), satisfies 

Kfifix) < Un(t,x) < a(t) + b{t)Ktip{x). 

Let u{t,x) := lim Un{t,x). Note that K^x) < u{t,x) < a{t) + b(t)Kt(p{x). 

Using the monotone convergence theorem, we have 

lim I dr I k(t-T, x,y)uf^(T,y)d}j.(y) = I dr I k{t-T,x,y)uP{T,y)d^i{y). 
"■^°°Jo Jm Jo Jm 

Since Un{t, x) satisfies 

u„+i{t,x) = Kt(p(x)+ f Krfix)dr+ f Kt^rKiT, 
Jo Jo 

we pass to the limit as n — > oo to obtain 



x)dT, (3.8) 



U{t, X) = Kfip(x) + I Krf{x)dT + I Kt-rU'\T, x)dT, 
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which shows that u{t, x) is a non-negative local solution of (ll.ll) - (ll.2l) for 
t e [0, To). 



Since a{t), b(t) are differentiable functions on [0, Tq), we see from (13.51) 
that for all t e [0, To), 

Mt, OIL < Mt) + b(t)K,(p\\oo < oo. 

The proof is complete. ■ 

Recall that, by Theorem 12. 3[ (ll.ll) - (ll.2l) does not have any essentially 
bounded global weak solution if a > /3 and p < and if f{x) ^ 0. How- 
ever, we can show that (|1.1I) - (|1.2I) possesses a essentially bounded global 
solution if p > for small functions / and (p (cf. [17] for Euclidean 
spaces). To do this, we need some integral estimates which are conse- 
quences of measure bounds for small and large balls. 

Recall that a measure on a metric measure space is upper a-regular if 
there exist some C, a > such that 

lx{B{x, r)) < Cr" (for all xeM,r> 0), (3.9) 

and is a-regular if there exists a constant C > such that for all .x; 6 M and 
all r > 0, 

C-^r" < iu(B{x, r)) < Cr" (for all xeM,r> 0). (3.10) 

It was shown in [llli Theorem 3.2] that if the heat kernel k satisfies (^5), (^6) 
with <l>2(i') satisfying 







s''-^(^2is)ds < oo, (3.11) 

then the measure p. is a-regular. Note that, by the monotonicity of 02, 
condition (13.1 II) implies that s"(^2(s) < C < oo for all s e [0, oo). 

Proposition 3.3. Assume that p is upper a-regular and xq is a reference 
point in M. IfQ<A\ < a and Ai -\- A2 > a, then there exists a constant 
Co > such that 

^ -My) < Co (for all x e M). (3.12) 



/ 

Jm 



d(y,x)^^[l d(y, xo)'^^] 



Proof. For each x e M, let Qi = {y e M : d(y, x) > d(y, xq)} and Q2 = 
M\Qi. Then 

X. d(y,x)Hl+d(y,Xo)'^]'^''^^^ " X d(y,Xo)Hl + d{y,Xo)'^'^^^^ 



and 



/ 



1 r 1 

-dp(y) < — . TT^dpiy). 



Ja, d(y,x)^^[l-^d(y,xo)^'^] Jm d(y,x)-^^[l d(y,x)'^'~] 

Routine estimates using upper regularity (13.91 ) now give uniform bounds on 
these integrals near xo and x (since Ai < a) and for large d(y, xo) and d(y, x) 
(since Ai -\- A2 > a), to give (13.121) . ■ 
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Proposition 3.4. Assume that fj. is upper a-regular and Xq is a reference 
point in M. IfQ<A\ < a and A2 > a, then there exists a constant Ci > 
such that 



j 

Jm 



^ ^My) • (3.13) 



Jm d(y, xY^ [1 + d(y, XqY^] 1 + dix, XqY' 

Proof. Fix X e M. If d(x,X(j) < 1, then (13.131) directly follows from 



(13.121) . since 



1 2Co 
-dfiiy) < Co < 



Jm d{y,xY'U +d{y,XoY'r'^^' l+d{x,XoY' 
Assume that d(x, xq) > I. If d(y, xq) > d(x, xo)/2, we have that 

1 C 

< 



\+d(y, XqY^ U+d(y, x^Y^-'^' ][\+d(x, XqY' ] 
where C is independent of xo,y. Using Proposition 13. 3 [ it follows that 

L{y,xo)>d{x,xo)/2 d(y,xY'[l +d(y,xoY']'^^^^^ 



c r 1 

1 + d{x, XqY^ Jd(y,xo)>d(x,xo)/2 d(y,xY'[l +d(y,xoY^-^']' 



< r- (3.14) 

1 + d(x, xoY' 

If d(y, Xq) < d(x, xo)/2, then 

2^1+1 



diy^xy"' < [d(x,Xo)-diy,Xo)y'' < [dix, Xo)/2]-'' < 



1 + d{x, XqY^ ' 

(3.15) 

and hence, 

r 1 C2 

Jd(y,xo)<d(x,xo)/2 d(y, xY'[l + d(y, xqY'] ~ 1 + dix, xqY' ' 
where we have used that f , „ ^ , , du.(y) < 00 as > a. 

JM [+d(y,xo) 2 

Adding (|3.14l) and (13.161) we see that (|3.13l) also holds if d(x, jcq) > 1 . ■ 
We now show the global existence of weak solutions for small cp and /. 

Theorem 3.5 (Global-existence). Let a > 13 > and suppose that the heat 
kernel k satisfies {k5), {k6) and that O2 satisfies l\3.11D . Let A> a and let xq 
be a reference point in M. The for each p > a /(a - j3) there exists 6 > 
such that if 

0<^{x),f{x)<-—^ 

1 + d{X, XqY 

then ( li.iD -( [772]) has an essentially bounded global solution. 
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Proof. Recall that conditions (^5), (^6) and (13.111) imply that jj. is a- 



regular. Let the map T be defined as in (13.61) : 

Tu{t, X) = K,(p{x) + I Krf{x)dT + I Kt-rU''(T, x)dT. 

Jo Jo 

For 6 > 0, let 5 e be the complete subset of the Banach space L°° ([0, oo) x M) 
given by 

ueV" ([0, oo) X M) : < u{t,x) < ^ J ^ . 

We will use the contraction principle to show that, for appropriately small e 
and 5, there exists a global solution in S^. 

For /i > a, we claim that there exists C2 > such that, for all < g{x) < 
6 /{I + d(x, xoY), we have 

c s 

Ktg{x) < for all X e M and all t > 0. (3.17) 

1 + d(x, Xq) 

To see this, let .x; e M. If d(x, Xq) < 1, then (13.171) is clear since 

Ktgix) = I k{t,x,y)g(y)diuiy) 
Jm 

- 1 1 , J -k{t,x,y)dn(y) 

Jm 1 + d(y, xqY 



< 6 I k(t, X, y)diu(y) < 6 

Jm 

25 

< -. 

1 + d{x, Xq)" 

So assume d{x, xq) > 1 . We have, using condition (^6) , 

r 6 

Ktg(x) < -k(t,x,y)diu(y) 

Jm 1 + d{y, XqY 

+ r T-77 -M,x,y)dfi(y)], (3.18) 

Ja^ 1 + d(y, XoY ) 

where Q.^ = {y e M : d(y, Xq) < d{x, Xo)/2} and Q.2 = M \ Q.i. For y e Q.u 
we have, noting from (13.111) that s"'(^2(s) is bounded, 

1 I d(y, x) ^ ^ 1 ^ d(y, x) )^' ^ ( d{y, x) 



t"IP "\ YIP I d(y,xY\ t'll^ I "\ t'lP 



C 2"C 
< < 



diy, xY d(x, XqY 

< , 

1 + dix, XoY 
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and hence, using that j^^^ < for A> a, 

r 1 1 (d(y,x)\ 2"+'c r 

c 
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< 



1 + d{x, xqY 



d(y, xoV 
(3.19) 



For y e Q.2, 



f L 



——k{t,x,y)dix{y) < 



< 



1 + d{x, xqY Ja.2 
C 

1 + d{x, XoY 



I k{t,x,y)dn{y) 

(3.20) 



using that A> a. Adding (13.191) and (13.201) . we see that (13.171) follows from 
(13.181) . proving our claim. 

Observe that by {k6) and (IXTTl) . 



I kiT,x, 
Jo 



y)dT < 



< 



< 



Jo r-IP H r'' 
/3 p 



.^) 



s"-f^-^^2(s)ds 



A — f s"~^-^(^2is)ds 
,xY-''Jo 



c 



d{y,xY~ 



(3.21) 



since 



Jo Jl 



S < +CX), 



using the monotonicity of O2 and (13.1 II) . 

Therefore, using (13.211 ) and (13.131 ) with A\ = a - /3 > and A2 = A > a. 



I Krf(x)dT = I I k(T,x,y)dT 
Jo Jm [Jo 

< 

Jm d{y,xY~'^(l+d(y,XoY) 



f(y)dM(y) 

6 



dfiiy) 



< 



C6 



1 + d{x, xoY~^ 



(3.22) 
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for all X G M and t > 0. Similarly, for m e 5e, we have that, using (13.131) 
with Ai = a - /3, A2 = p(a - yS) > or, 

r Kt-ruP(T,x)dT < f f k(t-T,x,y) ^ — —T-dfi(y)dT 
Jo Jo Jm (1 +d(y,xoy P)P 

C C sP 
~ Jm d(y, xT-P (1 + d{y, x,,r-P)p'^^^'' 

- ^''jMd{y,xY-n+diy,xor-P^p'^^^^ 

CgP 

- i+d(x,xor-'^ ^^'^^^ 

for all jc e M and ? > 0. It follows from (IXTTl) . (1X221) . (1X231) that if w g 5^, 
then 

^ , , C6 + CeP 

ru(t,x) < - + 



< 



< 



1 + d(x, XoT 1 + d(x, x^y-P 
Ci(6 + sP) 

1 + dix, xoY^I^ 
s 



1 + d{x, XqY-P 



provided that C\ {6 + sP) < s, in which case TS ^ c Sg. 

Next we show that T' is contractive on S g. Indeed, for m 1,^2 e S e, we 
have 

\Tuiit, x) - Tu2{t, x)\< \ I k{t - T, X, y) \u\{t, y) - u'^(t, y)\ diu(y)dT. 

Jo Jm 

Using the elementary inequality 

\aP - bP\ < pmsix{aP-\bP-^}\a-b\ for a,b >0,p > I, 
and the definition of S s, we obtain, using (13.211) and (13.121) . that 

\Tui(t, x) - Tu2(t, x)\ 



psP-^ 



Jo Jm [1 +d(y,xo)" > 



r c psp-' 

< C^ps'^'^WUl - U2\\ao- 

Thus if 6 is small enough to ensure that both C^psP~^ < 1 and CisP < s, and 
then 6 is chosen small enough so that Ci (6 + sP) < e, applying Banach's 



contraction principle to !F on the complete set S ^ implies that (II. 7|) and thus 
(il.ll) - (il.2l) has a global positive solution in Sg. ■ 
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4. Regularity 

In this section, we discuss the regularity of weak solutions. We show that 
weak solutions are Holder continuous in the spatial variable x if the source 
term / and initial value ip are both Holder continuous. We adapt the method 
used in 

In order to obtain the regularity of weak solutions, we need to assume 
that the function (D2 in condition (^6) satisfies the following assumption: 

Jr»oo 
s"^2{s)ds < 00, (4.1) 


where a is as in condition (^6). Since O2 is non-increasing on [0, 00), con- 
dition (|4~T]) implies that s^^"(^2{s) = o(l) as 5 ^ 00. 



Clearly, the Gauss-type function (D2 defined as in (12.41) satisfies condition 



4.11) for all y > whilst the Cauchy-type function O2 defined as in (12.5 



satisfies condition (14.11) for all y > 1 -I- a. 

Note that condition (14.11) is stronger than (13.111 ). and hence it implies that 
IX is a-regular. 

Proposition 4.1. Assume that ji is upper a-regular If ^2 satisfies d4.iD then, 
for all A e (0,1], 

£ d(x, y)'(:>2 (^^) diiiy) < C2t"'^'^ll' (for allxeM,t> 0) (4.2) 

for some constant €2- 

Proof. Let g(r) = r^(^2(jm) for r > 0. From (g^) 8(r) = o(r-") so, 
by a standard argument using a-regularity and integration by parts (see [0 
Proposition 4.1]), it follows that 

£d(x,yy^2{^^^yfi(y) = Jj{d(x,y))dii(y) 

<Ci r''\g'(r)\dr 
Jo 

= r y^'''^'- (^) ^ (^) '"'i 

As"'-'^'^^2is)ds+ I s"'-'^(-^2(s))ds 
^ Jo 



By an easy calculation, the last integral 

s"^'^ i-Q^'^is)) ds = - ^"+-^02(^)1^ + (a + ^) I s"'^-^-^Q)2(s)ds 
Jo 

Jr»oo 
s"^'^-^^2is)ds < C3 


using (14.11) . Therefore, 

^ j(x,j)''o(^j jMj) < C2t^"-'^'^ 
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as desired. ■ 

We now show the Holder continuity of weak solutions of (11.71) . 

Theorem 4.2 (Holder Continuity). Assume that (fi, f & L}{M) are Holder 
continuous with exponents 61,62 e (0, 1] respectively: for all Xi,X2 e M, 

\ip(x,)-ip(X2)\ < Csd{XuX2t, (4.3) 

1/(^1) -/fe)l < C(,d{xuX2t-, (4.4) 

where €5,0^ > 0. Assume that the heat kernel k satisfies (k5) - {kl) and 
that 02 satisfies ^4.71) with A = max {61, 62]. Let u(t,x) be a non-negative 
weak solution to ^i.il) - 4772]) that is bounded in (0, T) x M for some T > 0. 
Then u(t, x) is Holder continuous: for all Xi,X2 € M and all t e (0, T), 

\u(t, xi) - u(t, X2)\ < Cd{xi, X2f , (4.5) 

where 6 = 6icrl{6\ + vyS) and C > may depend on T but is independent of 
t, X. 

Proof. From (k6), (B31) . and (R^ . there exists C > such that for all 
t > and X e M, 

£k(t,x,y)\<p(y)-^ix)\dM(y) < Cjr'^/^ j^^^(^, 3;)^' O2 |^^jjM3^) 

< Ct^'"^. (4.6) 

By (11.71) it is enough to show that each of the functions Uq, U\,U2 is Holder 
continuous in (0, T) x M, where 

UQ{t, x) = K,(p{x), 

Uiit,x) = I Krfix)dT, 

Jo 

U2(t,x) = I Kt^rU'\T,x)dT. 



We first show the Holder continuity of uq. Indeed, for t > Q and X\,X2 € M, 
we see from {kl) that 



IMo(^-^i) - Uoit, X2)\ 



I {k{t, xuy)- k(t, X2, y))(piy)dp.(y) 
Jm 



< Lt ''d(xi,X2Tm\i 

< LMid{xuX2r-''° (4.7) 
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if f > d(xi,X2y°, where 5o > will be specified later on. On the other hand, 



if t < d(xi, X2y'', we have, using (^5), (1461) and (1431) . that 



\Uo(t,Xi) - Uo{t,X2)\ < 



I k{t, xi , y)i(p(y) - (fixi ))diu(y) 
Jm 

+ [(fixi) - (fiix2)] - I k(t, X2,y)i(p(y) - ip{x2))dn{y) 
Jm 



< ICf'^P + Csd{xuX2f 



< 



c[d{xuX2y''^'"^ + d{x^,X2f']. 



Combining this with (14.71) . it follows that 

|mo(?, -^i) - UoiU X2)\ < C [dixuX2r-''' + d{xuX2y°'''^ + d{xu X2f] 

< Cd{xuX2f'^'^'^'^\ (4.8) 

for all ? > and x\,X2 e M with d{x\,X2) < 1, where sq = cr/ (v + so 
that (T -vso = sqOi IP, and where we have used the fact that 6i > sqOi //3 for 
cr < 1 < vandyS > 1. 

Next we show the Holder continuity of Ui. As with (14.61) . we have from 
(k6), (E31) and (EJ) that 



Jm 



kiT,xuy) \fiy)-fixO\dM{y)<CT' 



which yields that, using (^5) and (|4.4I) . 



\Ui(t,Xi) - Ui{t,X2)\ 



k{T, xuy)(f(y) - f{Xi))dn{y) 



r [KJ{x,)-Krf{X2)]dT 

Jo 

Cdr f 

Jo Jm 

+t[fix,)-fix2)]- f dr [ k(T,X2,y)ifiy)-fix2))diu{y) 

Jo JM 

< 2C [ T^'"^dT + C(,td{xuX2f' 
Jo 



< 



Cf2//}+i _^Cetd(xuX2f' 

C [dixuX2y''''''^ + dixuX2y'^'''] 



(4.9) 



if t < d(xi,X2y\ where 5i > will be chosen later. 
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On the other hand, if t > J(xi,X2)*', and setting ti = d(xi,X2y\ we 
obtain, using (kl), that 

r [Krf(x,)-Krfix2)]dT < f cIt f \k{T, x^y) ' Kt, X2,y)\\fiy)\ dfi{y) 

Jtt Jti J M 



< f Lr-'dixuXjrWfWidr 



V- 1 



-d{xi,x2r\\f\\i 



< 



V- 1 

It follows from (14.101) and (1431) that 



d{^x\ , X2) 



(4.10) 



\Ui{t,Xi)-Ui{t,X2)\ < 



Jo 



f{x,)-KJ{X2)dT 



+ 



f{x,)-KJ{X2)dT 



< C [d{XuX2y'*'''^-"^ + J(JC1, X2)-^'^^^ + X2)-^'('"'^^^] 

< Cd{xuX2r^''^^^'^''^''''^ (4.11) 

if <i(jci, .x:2) < 1, where si = 0-/3/(62 + VjS) so that si + 8^62/ /i = 51(1 - v) + cr, 
and where we have used the fact that 

5l + 6*2 > 5i + 8162/ /3 

for cr < 1 < vandjS > 1. 

Finally, we show the Holder continuity of U2. Since u{t, x) is bounded on 
(0, T) X M, we see that 



JdT I k{t - T, X, j)m^(t, y)diu(y) < Crj. 
t-ij J M 



Hence, using {kl), we obtain 

\U2{t,Xi)-U2{t,X2)\ = 



dr I k{t - T,xi, y)u''{T, y)djjL{y) 



-fdrf 

Jt-ri JM 



k{t - T, X2, y)uP{r, y)d^(y) 

-J] Jm 

+ f dr f {k{t-T,xuy)-k{t-T,X2,y))u''{T,y)dniy) 
Jo Jm 

< 2C7] + L f dr f \t-Trdixi,X2TuP{T,y)dfi{y) 

Jo JM 



< C{r] + ri'~^d{xuX2r). 
Taking rj = d(xi, X2)°'^^, we thus have 

\U2(t, Xi) - U2(t, X2)\ < Cd{X\,X2T^^ . 

Combining (HTSl) . (14.111) and (14.121) . we conclude that 

\u{t,xi) - u{t,X2)\ < Cd{XuX2f"^'^^'^''^\ 



(4.12) 
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for all t e (0, T) and xi,X2 & M with d(xi,X2) < 1, for some C > 0, where 
we have used that 

010-/(^1 + v/3) < Cr/v < Cr(02 + /3)/(e2 + v/3). 

The proof is complete. ■ 

Finally, one may show that if the heat kemel k satisfies {k5), if \\f\\oo < °° 
and if (p{x) satisfies 

\Kt+s(p(x) - Kr(p(x)\ < C6 (for a\\t>0,xe M), 

then the essentially bounded weak solution u of (|1.7I) is Lipschitz continu- 
ous in time t on (0, T) x M, that is, 

\u{t + 6,x)- u{t, x)\ < Ci6 (t e (0, T),6>0,xe M). 

We omit the details, which are similar to the special case considered in [|9l. 

We note that, unlike the blow-up and the existence, the regularity of so- 
lutions is not related to the Hausdorff dimension a and the walk dimension 

References 

[1] C. Bandle, H. A. Levine and Qi S. Zhang, Critical exponents of Fujita type for 
inhomogeneous parabolic equations and systems, J. Math. Anal. Appl. 251, 624- 
648(2000). 

[2] M.T. Barlow, Diffusions on Fractals, Lect. Notes Math. 1690, Springer, 1998, 1-121. 
[3] M.T. Barlow and R.F. Bass, Brownian motion and harmonic analysis on Sierpmski 

carpets, Canad. J. Math. (4) 51 (1999), 673-744. 
[4] M.T. Barlow, R.F. Bass, Z.-Q. Chen and M. Kassmann, Non-local Dirichlet forms 

and symmetric jump processes. Trans. Amer Math. Soc. 361 (2009), 1963-1999. 
[5] M.T. Barlow and E.A. Perkins, Brownian motion on the Sierpmski gasket, Probab. 

Theory. Related Fields 79 (1988), 543-623. 
[6] Z.-Q. Chen and T. Kumagai, Heat kemel estimates for jump processes of mixed types 

on metric measure spaces, Probab. Theory Relat. Fields, 140 (2008), 277-317. 
[7] K. Dalrymple, R.S. Strichartz and J. P. Vinson, Fractal differential equations on the 

Sierpmski gasket, J. Fourier Anal. Appl. 5 (1999), 203-284. 
[8] K.J. Falconer, Semilinear PDEs on self-similar fractals. Comm. Math. Phys. 206 

(1999), 235-245. 

[9] K.J. Falconer and J. Hu, Nonlinear diffusion equations on unbounded fractal domains, 

J. Math. Anal. Appl. 256 (2001), 606-624. 
[10] H. Fujita, On the blowing up of solutions of the Cauchy problem for u, - A + u^^", J. 

Fac. Sci. Univ. Tokyo Sect. 1 13, 109-124(1966). 
[11] A. Grigor'yan, J. Hu and Ka-Sing Lau, Heat kernels on metric measure spaces and 

an application to semilinear elliptic equations. Trans. Amer Math. Soc. 355, 2065- 

2095(2003). 

[12] B.M. Hambly and T. Kumagai, Transition density estimates for diffusion processes 
on post critically finite self-similar fractals, Proc. London Math. Soc. (3)79 (1999), 
431-458. 

[13] J. Kigami, Analysis on Fractals, Cambridge University Press, 2001. 
[14] R.S. Strichartz, Differential Equations on Fractals : a Tutorial, Princeton University 
Press, 2006. 

[15] FB. Weissler, Local existence and nonexistence for semilinear parabolic equations in 
LP, Indiana Univ. Math. J. 29 (1980), 79-102. 



24 



FALCONER, HU, AND SUN 



[16] F.B. Weissler, Existence and non-existence of global solutions for a semilinear heat 

equation, Israel J. Math. 38 (1981), 29-40. 
[17] Q.S. Zhang, A new critical phenomenon for semiUnear parabolic problems, J. Math. 

Anal. Appl. 219(1998), 125-139. 

Mathematical Instttute, UNiVERsrrY of St Andrews, North Haugh, St Andrews, Fife 
KY16 9SS, UK. 

E-mail address: kj f@st-andrews .ac.uk 

Department of Mathematical Sciences, Tsinghua University, Beding 100084, China. 
E-mail address: huj iaxin@mail . tsinghua . edu . cn 

Department of Mathematical Sciences, Tsinghua University, Beding 100084, China. 
E-mail address: sunyh88©mails . tsinghua . edu . cn 



